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A spin nematic is a state which breaks spin SU(2) symmetry while preserving translational and 
time reversal symmetries. Spin nematic order can arise naturally from charge fluctuations of a 
spin stripe state. Focusing on the possible existence of such a state in strongly correlated electron 
systems, we build a nematic wave function starting from at — J type model. The nematic is a spin- 
two operator, and therefore does not couple directly to neutrons. However, we show that neutron 
scattering and Knight shift experiments can detect the spin anisotropy of electrons moving in a 
nematic background. We find the mean field phase diagram for the nematic taking into account 
spin-orbit effects. 
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I. INTRODUCTION 

When ordered, classical spin systems can arrange 
in a number of patterns, including (anti)ferromagnetic, 
canted, and helical structures. In addition to these, quan- 
tum mechanics allows the formation of a wealth of mag- 
netic phases for quantum spins not available to their 
classical counterparts. Due to its quantum numbers, 
detection of such order is often difficult: For instance, 
Nayak has considered a generalization of spin density 
wave (SDW) order, in which spin triplet particle-hole 
pairs of non-zero angular momentum condense with a 
modulated density 0. These states are characterized by 
spin currents rather than spin densities: thus, they do 
not couple at linear order to probes such as photons, 
neutrons, or nuclear spins. Only at second order do these 
phases couple to conventional probes, e.g. in two-magnon 
Raman scattering. Despite the challenges involved in 
their detection, subtle forms of magnetic ordering such 
as these may be necessary to explain phenomena such as 
the specific heat anomaly in the heavy-fermion compound 
URu2Si2@]> and the pseudogap regime in the cuprates Q. 

Promising materials for the observation of exotic mag- 
netic phases include systems with strong antiferromag- 
netic fluctuations such as the heavy-fermion compounds, 
the organic superconductors, and the cuprates. For in- 
stance, the cuprates in the absence of carrier doping are 
antiferromagnetic Mott insulators at low temperatures. 
As carriers are introduced through doping, the nature 
of the magnetic order evolves until, for optimally doped 
and overdoped samples, the system becomes a metallic 
paramagnet. In between these two limits, the under- 
doped cuprates have been argued to have spin glassQ 
and stripe phases 0. The proximity between Mott insu- 
lator and superconducting phases in the cuprates makes 
them ideal systems to study the hierarchy by which the 
broken symmetry of Mott insulators is restored [|| QJ • 

In this paper we will explore the possibility of detec- 
tion of spin nematic order, a different quantum magnetic 
phase, in strongly correlated electron systems @, U Ne- 
matic order has been proposed as a state originating from 
charge fluctuations of stripe order^,^]. A spin stripe 
is a unidirectional collinear spin density wave (SDW), 



and it consists of antiferromagnetically ordered domains 
separated by anti-phase domain walls, across which the 
direction of the staggered magnetization flips sign. The 
order parameter is 



S(r) = *e 



iK a r 



-iK,T 



with K s = (tt, 7r) + 5 corresponding to antiferromagneti- 
cally ordered stripes with period 27r/|<5|. Here, the com- 
plex vector $ = e %6 n takes its value within the manifold 
of ground-states S\ x S2/Z2; the Z2 quotient is necessary 
not to overcount physical configurations, as the transfor- 
mation e 10 — ► —e l6 , n — * — n does not modify 3?. The 
real vector n gives the direction of the staggered mag- 
netization in the middle of a domain, while the phase 
factor e w specifies the location of the domain walls. A 
shift in 9 by 2tt translates the system by the periodicity 
of the SDW, and thus leaves the system invariant. As- 
sociated with collinear SDW order is charge order due 
to the modulations in the amplitude of the local spin 
magnetization [T^ . which can be described by a general- 
ized CDW order parameter 



Sp(r) = <pe 



2iK s 



-2iK s 



for some SU(2) invariant observable p, not necessarily the 
electron density. In the stripe picture, this CDW usually 
arises from the accumulation of holes at the domain walls. 

Stripes were first observed in elastic neutron scat- 
tering experiments on the spin-1 nickelate insula- 
tor La2_ 3; Sr 2 ;Ni04, and coexistence of stripes with 
superconductivity was first observed in underdoped 
Lai «_»Ndn /i.Sr^CuO/i.[l3|. where the unidirectional 
character was demonstrated by transport^] and 
photoemission [l5j measurements. In YBa2Cu30g.35 and 
Lai.6-xNdo.4Sr 2; Cu04, it is observed that SDW order 
is first destroyed by spin fluctuations^}. In this case, 
memory of the charge modulation can remain, even af- 
ter averaging over the spin direction, resulting in a spin- 
invariant CDW phase, whose presence may explain recent 
STM measurements in Bi 2 Sr 2 CaCu20 8+ 5 [l7llTs|. 

For other materials, however, or in other regions of the 
phase diagram, symmetry may be restored in a different 
order. In particular, Zaanen and Nussinov [lflL Hll | pro- 
posed that many experimental features of LazSri-xCuC^ 



2 




FIG. 1: Charge fluctuations of a spin stripe, provided that 
the domain walls (curves) maintain their integrity. A spin 
nematic state is a linear superposition of fluctuating do- 
main configurations such as these. Within each domain, the 
staggered magnetization (arrows) is well-defined, and it flips 
sign across the anti-phase domain walls. Due to fluctua- 
tions, translational symmetry is restored to the system, and 
the magnetization has expectation value zero at every point. 
However, SU(2) symmetry is not restored, as there is a pre- 
ferred axis in which spins align, modulo a sign. [TH ITU. 



can be explained by assuming that the spin stripe order 
is destroyed by charge fluctuations. In this picture, dy- 
namical oscillations in the anti-phase domain-walls of a 
spin stripe lead to a restoration of translational symme- 
try and a loss of Neel order. However, although both 
charge and spin order seem to be destroyed in this pro- 
cess, 5p = and S = 0, full spin symmetry need not 
be restored: So long as neither dislocations nor topo- 
logical excitations of the spin proliferate, the integrity 
of the domain walls allows the staggered magnetization 
on each oscillating domain to be well-defined. Thus, al- 
though the local magnetization does not have an expec- 
tation value, the magnetization modulo an overall sign 
does. This is nematic order, in which S and — S are iden- 
tified, see Fig. ^ The order parameter can be chosen to 
be (S a SP - S 2 6 al3 /3) oc (n a n - l/in 2 5 al3 ) j= 0, which 
has a non-zero expectation value. Because translational 
invariance is restored in this process, the nematic order 
parameter is spatially uniform, instead of being modu- 
lated by some multiple of the SDW wave vector. In ad- 
dition, we expect the nematic to be uniaxial, with a sin- 
gle preferred axis n (mod Z 2 ) inherited from the nearby 
collinear SDW. 

Direct observation of spin nematic order through con- 
ventional probes is difficult. For instance, neutrons do 
not couple to nematic order, which is a spin-two opera- 
tor. Similarly, nematic order is translationally invariant, 
and does not give Bragg peaks in X-ray experiments. In 
principle, two-magnon Raman scattering can probe ne- 
matic order, but in practice it is difficult to separate the 
contribution due to the nematic from the creation of two 



magnons[10j. Hence, although spin nematic order can 
have important experimental consequences, e.g. for anti- 
ferromagnetic correlations in magnetic field experiments 
on superconducting samples [Io|. its direct detection re- 
mains a challenge. 

The existence of other stripe liquid phases, different 
from the spin nematic treated in this paper, as well as 
proposals for their detection, are discussed in Ref. 0. In 
particular, the nematic state described there originates 
from fluctuations of a unidirectional CDW that restore 
translational invariance but, by maintaining a memory of 
the original orientation of the CDW, break the rotational 
symmetry (point group) of the lattice. Hence, unlike 
the spin nematic, this "charge nematic" is SU(2) spin 
invariant, and only breaks a discrete group. 

II. SPIN NEMATIC ORDER PARAMETER 

A spin nematic is a state that breaks spin SU(2) sym- 
metry without breaking time reversal invariance |g, |20j . 
The presence of spin nematic order can be observed in 
the equal time spin-spin correlator 

(S a ( ri )^(r 2 )) = C(r 1 ,r 2 )^ + e Q/5 ^(r 1 ,r 2 ) 

+Q aP (r 1 ,r 2 ). (I) 

This expression corresponds to the SU(2) decomposition 
(I) (1) - (0) sym 8 (l)asym © (2)s ym . We consider the 
three terms appearing in eq. in turn. The scalar func- 
tion C is explicitly spin invariant. It contains important 
information regarding charge order, but it does not help 
us in defining a nematic phase. On the other hand, the 
pseudovector function A = (Si x S 2 ) gives a measure 
of the non-collinearity of the spin vector field. However, 
in the case at hand where the nematic state originates 
from charge fluctuations of a collinear SDW phase, we 
expect A to vanish. This is supported by the fact that, 
from the point of view of the Ginzburg-Landau free en- 
ergy, the pseudovector A cannot couple linearly to any 
function of the SDW order parameter <&. As an aside, 
we note that for systems with spin exchange anisotropy 
of the Dzyaloshinskii-Moriya (DM) form, the DM vec- 
tor will couple linearly to A, as expected from the weak 
non-collinearity (canting) in such systems. However, the 
expectation value of A in this case comes from explicit 
breaking of the spin symmetry. 

Thus, all information of interest to us is contained 
in the symmetric spin-2 tensor Q a @ . We define a sym- 
metrized traceless spin correlator Q a/3 , 

g^(n,r 2 ) = ±(§°§i + S?S«) - ^Si ■ S 2 (2) 

whose expectation value yields Q a/3 directly, 

Q^(r 1 ,r 2 ) = (Q^(r 1 ,r 2 )). (3) 

Starting from a SDW state, as domain wall fluctua- 
tions grow to destroy charge order, translational invari- 
ance is restored to the system, insuring that Q a/3 (ri,r 2 ) 
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is independent of the center-of-mass coordinate, i.e. 
Q a P(ri, r 2 ) = Q a ^(r 1 + R, r 2 + R) for any displacement 
R. This fact alone signals the breaking of spin sym- 
metry, since the choice of a non-trivial (i.e. not pro- 
portional to S ab ) tensor Q a ° has been made across the 
system. Thus, while Q a/3 (ri,r 2 ) decays exponentially 
with distance |ri — r 2 | due to the absence of long-range 
Neel order, the onset of nematic order is reflected in the 
translationally-invariant expectation value in the matrix- 
valued function J3J). Since the original SDW state has a 
single preferred spin direction n, the ensuing nematic or- 
der will be uniaxial, 

Q^(r 1 ,r 2 ) = f(v 1 -r 2 )Qf, 

(4) 

Qf=(n a nP-5 a < , /3)S. 

Here, we have decomposed the order parameter into three 
component objects: a function / describing the internal 
structure of the nematic; the unit vector director field n; 
and the scalar magnitude S. As seen explicitly from |J3J, 
the function / is parity-symmetric, /(r) = /(— r). As 
will be shown below, / is dominated by the short range 
antiferromagnetic correlations between spins, and there- 
fore has a large contribution at the wave vector (7r,7r). 
By definition, we choose this contribution to be positive, 
/(■7T, 7r) > 0. With this convention, S > corresponds 
to a Neel vector that is locally aligned or anti-aligned 
with the director field n (sometimes referred to as the 
N + phase in the literature of classical liquid crystals, see 
e.g. Ref. I21T) . while S < corresponds to a Neel vec- 
tor that is predominantly perpendicular to n (the N - 
phase). We will show that, at low temperatures, S > 
due to the local antiferromagnetic correlations whereas, 
for anisotropic systems at high temperatures, a phase 
with S < is possible. Finally, we note that the uni- 
directional SDW state (stripe phase) breaks the discrete 
rotational symmetry of a tetragonal lattice. This sym- 
metry may be restored when charge fluctuations destroy 
the SDW state to form the nematic. In this chase /(r) 
will be symmetric under rotations on the plane by 7r/2, 
r — * R^/2 V - However, if a memory of the orientation 
of the stripes survives the domain wall fluctuations, / 
will not have such symmetry, yielding a "nematic spin- 
nematic" , i.e. a translationally invariant system that is 
anisotropic in real space and in spin space. 

Sections IIIII and |IVI will be devoted to understanding 
the behavior of the three component fields of the nematic: 
/, n, and S. Then, in Section ]V\ we will explore how this 
detailed knowledge can be used in experimental searches 
for nematic order. 



III. NEMATIC WAVE FUNCTION 

In order to explore the possible symmetry properties of 
the function /, we study its short wave length structure 
by explicit construction of a nematic operator on a small 
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FIG. 2: Spectrum of t-J model on a 2 x 2 plaquette. Adapted 
from Q. 

cluster. As is well known, it is impossible to describe 
nematic order in terms of a single spin-1/2 particle, as 
the identification of "up" and "down" results in a trivial 
Hilbert space for the spin degree of freedom. Another 
way to see this is that a spin-2 operator has a vanish- 
ing expectation value with respect to any spin-1/2 state, 
resulting in the identity Qf- = whenever i = j. 

This limitation can be overcome by coarse-graining a 
group of spins and constructing a nematic wave func- 
tion out of them. In order to preserve the underlying 
rotational symmetry of the system, we carry this out on 
square 2x2 plaquettcs of spins, and use energetic con- 
siderations to find the states most likely to contribute 
to magnetic order. We take the t-J model as a starting 
point, analyzing the low energy Hilbert space in a man- 
ner similar to the projected SO (5) approach of Zhang et 
al. [2^ or the CORE approach of Altman and Auerbach 
[23| . In these analyses, the lattice is first split into pla- 
quettes. On each plaquette, the Hamiltonian is diagonal- 
ized, and the m lowest energy states are kept, where 
v € {1, . . . , m} and i labels the plaquette. For a lattice 
composed of N plaquettes, this allows one to define a pro- 
jected subspace M of the Hilbert space spanned by m N 
factorizable states of the form |VVl)i ® |Vv 2 )2 ■ • • \<Pv n )n- 
We assume that the ground state is well-contained in M. . 

Figure [3 reproduces results in [23| for a t-J model on 
a 2 x 2 plaquette. The lowest energy bosonic states are, 
at half filling, the S = ground state |f2) and the S = 1 
magnon triplet ijjfi), and for a plaquette with two holes, 
an S — state 6'|f2). In addition, there are two low- 
lying S — 1/2 fermion doublets with one hole; however, 
unbound holes are dynamically suppressed, as supported 
by DMRG calculations on larger lattices, and we shall 
exclude the one-hole sector in what follows. Then, to 
lowest order in the analysis, we only keep the low- lying 
bosonic states, in terms of which a general low-energy 
plaquette state can be written as 

\il>)i= (e + m^+dbfjpU. (5) 

It is useful to get some intuition regarding the wave 
functions JSJ. Introducing the total spin and staggered 
spin operators on a plaquette, S = Si + S 2 + S3 + S4 and 
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N = Si — S2 + S3 — S4, as well as the DM-type vector, 
D = Si x S 2 - S 2 x S 3 + S 3 x S 4 - S 4 x Si, we find 

■^V a ^ H~ ; 



A, 



*(^a 



up to positive multiplicative factors and up to terms lying 
outside of the projected low-energy space. In particular, 
we find the expectation values on a single plaquette 

(Si) ~ m* x m, 
(Nj) ~ Re(s*m), 
(6*) ~ Im(s*m), 
(bi) ~ s*c, 
(N,&j) ~ m*c. (6) 

The last two expressions describe local singlet and triplet 
superconductivity, respectively. 

One can build a nematic on a cluster of plaquettes by 
choosing a factorizable state, in which s = and m is a 
constant real vector on every plaquette, 



\^ct) = U(m a il l + cb\)\n) t 



(7) 



The constraint that m is real insures that no long range 
ferromagnetic or Neel orders develop, whereas nematic 
order does due to the SU(2)-symmetry breaking choice 
of the vector m. Note from © that IV'fact) is also 
a triplet superconducting state (except at half-filling, 
where c = 0). Order of this type is found, for instance, 
in the triplet superconducting state of quasi-one dimen- 
sional Bechgaard salts, where the triplet superconduct- 
ing order parameter is constant along the Fermi surface 
due to the splitting of the Fermi surface into two disjoint 
Fermi sheets j^. 

On the other hand, in applications to materials such 
as the high T c cuprates, we would like to introduce a 
nematic state that is a singlet superconductor, instead 
of triplet. For this, the use of non-factorizable states is 
necessary. For instance, introducing a local angle variable 
9i £ [0, 2tt) on each plaquette, consider the state 



J(d6 1 

IB- 



d9o 



(8) 



; cos(Q 



where m is a constant real vector, and Q is the wave 
vector of the underlying SDW. If the angle Qi were held 
constant across the lattice, long range SDW order would 
ensue. In contrast, by integrating independently over 
the 6i at different sites, we introduce charge fluctuations 
that average out the local magnetization to zero. Hence, 
the state JSJ has restored translational and time-reversal 



symmetry, with only singlet superconductivity and ne- 
matic order surviving. We note that © ignores correla- 
tions between spin degrees of freedom, controlled by v a , 
and charge degrees of freedom, controlled by w . More 
complex nematic wave functions that take this effect into 
account are given in Appendix 1X1 On the other hand, in 
practical calculations, one may consider a slightly sim- 
pler wave function than JSJ by replacing the 9i by local 
Ising variables a i — ±1 on each plaquette. This leads to 
the wave function 



IV>2> 



E 



{CT1,CJ2...} — ±1 



n 



{8 + a i m a P ai + cb})\{l) i 



(9) 



As before, this state has time-reversal and translational 
symmetries restored, with only spin nematic and super- 
conducting orders surviving. Finally, note that in order 
to produce a nematic state without any type of super- 
conductivity (singlet or triplet), it is necessary to intro- 
duce another fluctuating Ising variable which multiplies 
the term cb\ in (JSJ, thus randomizing the relative phase 
between all three components of the wave function. 

Despite the fact that we have considered many different 
wave functions, depending on the possible coexistence of 
nematic order with different types of superconductivity, 
the analysis above allows us to make strong predictions 
on the dominant short wave length dependence of /. This 
is because, of all the low energy plaquette states kept in 
the projected Hilbert space, only the triplet tjjft) breaks 
SU(2) symmetry. Hence, only this state can contribute 
directly to the nematic order parameter. There are six 
different pairs i ^ j of sites on a 2 x 2 plaquette, and the 
most general spin-2 operator P a/3 on a plaquette can be 
written as a linear combination of the six "link" operators 



P 



a/3 



(10) 



It is useful to introduce an inner product for real func- 
tions ctij on the links of a plaquette, according to 
(a,/3) = X3rjj.ji.gQ otijflij. With this, we can write a 

normalized basis of plaquette operators P v , correspond- 
ing to 6 linearly independent functions which satisfy 
(a v ,a") = S ,,l/ . We choose a basis of eigenoperators of 
the symmetries of the plaquette, composed of operators 
with p x , p y , d x 2 +y 2 and d xy symmetries, see Figure|21 in 
addition to two operators with s-wave symmetry, 



pa/3 



1 [Q?i 



Q23 



Q41 



(11) 



It is easy to see that among the six operators ffTTTfl . 
only those with s-wave symmetry get a non-vanishing 
expectation value with respect to the state JSJ. What's 
more, the linear combination P SF = (P Sl + >/2P S2 ) / y/3 
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FIG. 3: Basis of plaquette nematic operators with well- 
defined symmetry properties. The solid lines denote links 
with weight +1, dashed lines have weight — 1. Thus, for in- 
stance, P S1 = (Q12 + Q23 + Q34 + Q41) /2, which in terms of the 
Q.ij vectors reads a si — (ei 2 + e 2 3 + £34 + e4i)/2. The other 
basis vectors are a S2 — (en+^i)/^, oP x — (623 — e4i)/-v/2, 
a Py = (e i2 - e 34 )/\/2, a d * 2 +« 2 = (e i2 - e 2 3 + e 34 - e 4 i)/2, 
and a*"* = (eia - e 24 )/\/2. 



has a vanishing expectation value, so that it is possible 
to write an orthonormal basis of operators in which only 
the basis operator, 



P SA = ^(-V2P S1 +P S2 ), 



(12) 



has a non-zero expectation value. This can be understood 
as a consequence of local antiferromagnetic correlations, 
since can be expressed in terms of the plaquette 
staggered magnetization N as 



1 / i) a/3 
PfP = -L= [N a N^-N 2 — 
SA 2V6 \ 3 



Using the relation 



n 



(13) 



(14) 



which holds for i,j on the same plaquette, we see that the 
dominant short-range contribution to the nematic order 
parameter is 



(15) 



This is of the form J3J with 



/(k) = a SA (k) 
4 

:(cos k x cos k y — cos k x — cos k y ), (16) 



.a/3 
10 



Note that we can relate the real vector m in Eqs. © 
and to the director field n of the nematic, appearing 
in through 



(P°f ) cx m a mP - m 2 S al3 /3, 



(17) 



from which we conclude that nam. This can be used 
to constrain the sign of S, 

S = |nV<i^f) > 0. 



N : S<0 




FIG. 4: Mean field phase diagram for a nematic without 
anisotropy in the spin exchange interaction. The ordered 
phases N ± axe uniaxial. All transitions are first order. 



Therefore, at low temperatures, where the state of the 
system is dominated by the low energy plaquette states, 
S is positive. 



IV. MEAN FIELD ANALYSIS 

We now study the finite temperature phase diagram 
of the spin nematic. We assume that superconductivity 
is either present as a background phase throughout the 
entire region of the phase diagram that we study, or not 
present at all. Hence, we do not include explicitly the 
interaction between superconducting and nematic order 
parameters. As is well known (see, e.g. Ref l2o|) . symme- 
try allows the inclusion of cubic terms into the Ginzburg- 
Landau (GL) free energy of a nematic order parameter. 
Thus, the most general GL free energy for a spin-isotropic 
system is, up to quartic order, 



^Q.iso = pTrQl ~ lTrQ 3 + 6TrQ 4 . 



(18) 



By the identity (TrQg) 2 = 2TrQ^, which holds for any 
3x3 traceless symmetric matrix, we do not include the 
term (TrQg) 2 in (JTHJ}. The P hase dia 

gram is shown in 

Fig. 0] At high temperatures, (5 is large and positive, 
and the system is disordered, Qo = 0. On the other hand, 
as temperature is reduced and the value of (3 decreases, 
there is spontaneous breaking of spin SU(2) symmetry 
as the system undergoes a first order transition into the 
nematic state. The nematic in this case is uniaxial as 
desired (see sectionsllland lll|) : in order to stabilize biaxial 
nematic order, terms of order Qq and Qq would have to be 
added to the free energyj2f|. We note that, as discussed 
in sectionslllland lTTll S > in the low temperature phase. 
Thus, the cubic coefficient 7 must be positive. 

How is the above analysis modified by the presence of 
spin-orbit effects? For definiteness, we concentrate on 
the case of the cuprate La x Sri_ 2; Cu04 in the low tem- 
perature orthorhombic (LTO) phase. This compound 
displays strong evidence for fluctuating stripe order in 
its underdoped regime |2fij. and the spin exchange con- 
stants J a/3 are known in detail from neutron scattering 
experiments in undoped Lao CuOa [27I |2S| . The analysis 



presented here can be easily generalized to other materi- 
als. In La 2; Sr 1 _ x Cu04, spin-orbit effects are small: The 
anisotropic part of the spin exchange interaction J Q/3 is 
less than 10~ 2 of the isotropic part in undoped La2Cu04 
[27L |28| . Yet, at low temperatures, this anisotropy leads 
to a preferred direction for the Neel order and to weak 
ferromagnetism |29j |. Thus, although the anisotropy is a 
low energy effect, playing a weak role on the onset and 
magnitude of the nematic order parameter, it may ulti- 
mately fix the preferred spin orientation for the nematic. 
This aspect of the interplay between J and Q will be 
especially important when we try to separate their con- 
tributions to the anisotropy in the spin susceptibility (see 
Eq. ||l>j! below). 

The most general form of the spin-exchange interaction 
for spins on nearest neighbor Cu sites is 

n = Y^-K'/si's; (19) 

(ij) 
(ij) 

where is a traceless symmetric tensor and D is the 
Dzyaloshinskii-Moriya (DM) vector. Neither Jo nOT Js 
depend on the bond (ij). It is convenient to work with 
the principal axes of the system: the vectors a and c lie 
on the CuC>2 planes at 45 degrees from nearest neigh- 
bor Cu-Cu bonds, and the vector b is normal to the 
CuC>2 planes, see Fig. |SJ In this basis, J s is diagonal, 
J 1 + J s = diag( J aQ , J bb , J cc ) with J bb < J aa « J cc H3- 
For simplicity we take J aa — J cc in what follows, cor- 
responding to easy plane antiferromagnetic interactions 
along the a-c plane. On the other hand, the DM vector 
Djj points along ±a, with the sign on each bond (ij) 
given by the staggered pattern shown in Fig. [SJ Thus, 

/ Jo + A \ 

Jij = J - 2A ±D , 
V tD Jo + A J 

where A = (J aa - J bb )/3 > 0. A useful measure of the 
relative importance of the anisotropies A and D is given 

n 2 

by the ratio x = -t~A' wri ich is approximately equal to 
one in the LTO phase of La 2 Cu0 4 HI Hi- 
lt is instructive to consider first the effects of spin 
anisotropy on an antiferromagnet. The GL free energy 
Fjy of a Neel order parameter N is, to lowest order in 
spin-orbit coupling, 

\i 

F N = -AN T J S N+— (D ■ N) 2 + /iN 2 + idST 4 . (20) 
Jo 

The term D • N is a pseudoscalar and is forbidden by 
parity and time-reversal symmetries. The first two terms 
include the effects of anisotropy, and are both quadratic 
in the spin-orbit coupling. The coefficient A is positive, 
as necessary to capture the tendency of spins to order 
along the easy plane at low temperatures. Similarly, the 
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FIG. 5: The Dzyaloshinskii-Moriya vectors (dashed arrows) 
and ground state spin orientation (solid arrows) on a single 
CUO2 plane in the LTO phase of La2Cu04. The black dots 
denote Cu sites. The DM vectors D^j point along ±a and 
are staggered between adjacent bonds. This, combined with 
easy plane anisotropy, yield a weakly canted AF order, where 
the spins on the Cu sites have a large staggered component in 
the ±c direction, and a small uniform component out of the 
plane (in the b direction, not shown). 

DM vector induces a weak canting for staggered spins 
that are perpendicular to it. This pushes N towards the 
plane normal to D, which requires A' > 0. Thus we see 
that J s chooses an easy plane for the spins, and the DM 
vector selects a preferred direction, c, within the easy 
plane [2^|, see Fig. EJ The ratio A'/A is unimportant in 
this case. 

We now turn our attention to the nematic. In addition 
to the usual GL free energy for a nematic order param- 
eter, eq. (|18fl. the explicit symmetry breaking due to 
anisotropy in the spin exchange can be taken into ac- 
count, to quadratic order in the spin-orbit interaction, 
by adding the terms, 

Fq — Fq >1so + Fc^anis 

F Q ,anis = -aTr(J s g ) + ^D T Q D. (21) 

Note that, unlike the Neel case l|2Uf) . the order parameter 
Q couples linearly to the anisotropy. Thus, the symme- 
try is broken explicitly, and strictly speaking there is no 
disordered phase with Qq = 0. This, however, does not 
preclude the existence of crossovers in the order param- 
eter, or even discontinuities at first order transitions, as 
we traverse the phase diagram. For weak anisotropy the 
order parameter will be extremely small at high tem- 
peratures, and very sharp crossovers will be observed. 
Another consequence of the linear coupling is that, un- 
like antiferromagnetic order which always lies on the easy 
plane of J Q ^, uniaxial nematic order can point either in 
the direction of maximal coupling of J a ^ , or in the direc- 
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FIG. 6: Mean field phase diagram including anisotropy in the 
spin exchange J Q/3 , for the case J bb < J aa = J cc , D = ±a. 
All transitions are first order. Phase N + is characterized by 
S > and n = c, phase N~ by S < 0. The director n in 
phase N~ depends on the parameter w = a'x/a: for w < 1, 
n = b, while for w > 1, n = a. The parameter w is material- 
dependent and is unlikely to change across the phase diagram. 
Thus, only a single type of N~ phase is accessible within a 
given material. The dashed arrow shows a possible trajectory 
where a first order transition is crossed as temperature is in- 
creased; alternatively, for smaller values of 7, the first order 
transition may be avoided and strong crossover behavior may 
be observed instead. 
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FIG. 7: Value of 5* along the dotted trajectory in Fig. |S| 
The parameters a A = 0.1, w < 1, and 7 = 5 = 1 are cho- 
sen for concreteness; /3 increases with temperature. As the 
anisotropy A is decreased, or for trajectories that start deep 
inside the N + phase due to a larger value of 7, the value of 
\S\ in the high temperature phase is reduced. Similarly, for 
small values of 7, the first order transition may be avoided 
and instead a strong crossover within the N~ phase may be 
observed. 



V. DETECTION 



tion of minimal coupling of J a @ , depending on the sign 
of S. For the state with nematic order and S > 0, we 
expect the director field n to lie on the easy plane and be 
orthogonal to D, i.e. be parallel to the direction that the 
SDW state would take in the absence of charge fluctua- 
tions. This constrains the linear coefficients appearing in 
(|21|l to be positive, a > and a' > 0. 



In order to obtain the mean field phase diagram for a 
spin nematic, we compare the minima of F for the direc- 
tor field n pointing along the various principal axes e,-. 
The result is shown in Fig.[S] The low temperature phase 
is a uniaxial nematic characterized by S > and 
n = c. The high temperature phase N~ is also uniaxial, 
but it has S < 0. The director field in N~ depends on the 
relative strength of the two anisotropy terms, w = a'x/a. 
This quantity is material dependent, and is unlikely to 
change significantly over the phase diagram (except, of 
course, across a structural phase transition), so that only 
one of the following scenarios should be observed within 
a given material: For w < 1, n = b, whereas for w > 1, 
n = a. The coefficient (3 increases with temperature, pos- 
sibly tuning a first order phase transition between phases 
N + and N~, or otherwise moving the system through a 
sharp crossover within the N - phase. In either case, the 
temperature dependence of the order parameter leads to 
strong experimental signatures discussed in Section 
Figure [7| shows the value of S as we move across the 
N + /N~ phase transition along the dashed line in Fig. 



In the presence of a spin nematic order parameter 
Q°j = Q a/3 (ri,Tj) of the form J3J|, symmetry allows the 
term 



Hi 



<ij °i S j 



(22) 



to enter the Hamiltonian, which can be thought of as a 
local spin-spin interaction mediated by the nematic or- 
der. The coefficient g is positive, as necessary for the sta- 
bility of the nematic order. When the system develops 
long-range nematic order, electrons move in a nematic 
background which acts like an effective anisotropic spin 
exchange. This leads to anisotropy in the spin response 
function. In the random phase approximation (RPA), 
when an external magnetic field B is applied to the sys- 
tem, H = Ho + Hint — Yl% • Sj, the electrons see an 
effective field 

B c Vk) = B*(k) + 5 tr 3 (k)(^(k)), (23) 

leading to the dynamic spin response 

X^ A (k,o;) =xo(k,w) [(i 2x2 +.gg(k) X o(k, w )) _1 ]t24) 

Spin nematic order thus induces anisotropy in the spin 
susceptibility. We introduce a tensor S"^ 3 , defined by 



Experimental measurements of from polarized neu- 
tron scattering experiments (for arbitrary wave vector k) 
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can be used to compute S a/3 , which is the natural object 
to consider when studying nematic order. In the current 
approximation, we find 



(0,71) 



S^ A (k, W ) = gQ^(k). 



(25) 



Similarly, Knight shift measurements give access to the 
local static spin susceptibility, i.e. the susceptibility in- 
tegrated over all wave vectors k. To linear order in the 
nematic order parameter Q, the anisotropy in the Knight 
shift over the various principal axes a is proportional to 



lhk [ Xo (k, W = 0)] 2 S^(k,c = 0), 



(2tt 



where no summation over a is implied. 

Even in the absence of nematic ordering, spin-orbit 
coupling leads to anisotropy in the antiferromagnetic ex- 
change Jff , which enters the Hamiltonian in a term of 
the form 122|) with gQ — ► — J, see eq. I|19[) . The complete 
expression for S Q ^ is therefore 



S^ A (k,«) = ff Q^(k)-J^(k). 



(26) 



Note that only the traceless symmetric part of Jjj con- 
tributes to 5. We face the challenge of untangling 
the contributions to i|2rj|) due to the anisotropy in J 
from those coming from the presence of nematic or- 
der. For this, the analysis of Section 11111 is pivotal, 
in particular eq. (|16|l . which gives the dominant wave 
vector dependence of the nematic contribution to the 
anisotropy in the spin response function, see Figure [5] 
This can be combined with detailed knowledge of the 
form of J s [13, |H, J s (k) = 7y(k)A diag(l, 1, -2). As 
before, A = (J aa - J bb )/3 > 0, and we assume that 
anisotropy is small beyond the nearest-neighbor range to 
set ?7(k) = 2(cosfc. E + cosk y ). Thus, the contributions 
due to Q and J can be distinguished by the wave vec- 
tor dependence of the signal, as one is proportional to 
/(k) = -^(cosk x cosfc y — cosk x — cosk y ) and the other 
to f](k) — 2(cosfc x + cosky). For example, experiments 
measuring the susceptibility at the points k = (tt, 0) and 
(0, 7r) are sensitive to Q, but not J, see Fig. |S] 

For completeness, we evaluate the expression (|26[) 
in the various phases shown in Fig. for underdoped 
La :E Sri_ 3 ;Cu04. In the low temperature phase N + , we 
expect, 

S aa = -g\S\f (k)/3-A n (k), 
E bb = - 5 |5|/(k)/3 + 2A„(k), 
S cc = 20|S|/(k)/3-At,(k), 

In the high temperature phase N - , we must consider two 
separate scenarios. If w > 1, then n = a, and 

= -2 fl |S|/(k)/3-Arj(k), 
E bb = g\S\f (k)/3 + 2 A V (k), 
S cc = g\S\f(k)/3-A V (k), 
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FIG. 8: Wave vector dependence of /(k) corresponding to the 
plaquette nematic operator P SA . This is expected to be the 
dominant contribution to /, and enters the anisotropy of the 
spin susceptibility through Eq. 1261 The "+" and "-" regions 
show the sign of /(k) on the Brillouin zone. 



whereas, if w < 1, then n = b, and 

~ cc = g\S\f(k)/3-Ar,(k), 

Thus, if the system goes through a phase transition be- 
tween phases N + and N - , the discontinuity in S, shown 
in Fig.0 together with the change in the director field n, 
would give a very clear experimental signature of nematic 
order in Knight shift and polarized neutron scattering ex- 
periments. 

Finally, we consider the prospect of observing direc- 
tor density waves (DDW), the Goldstone modes corre- 
sponding to the nematic state, in neutron scattering ex- 
periments. To do this, we study the pole structure of 
the RPA result Expanding J2H) about k = and 

u> = 0, we find two degenerate DDW modes, with speed 
«DDW oc y/g\S\, 



v 2 k 2 ' 
y DDW K 



(27) 



with _L labelling either of the directions perpendicular 
to the director of the nematic. Equation l|27|l is con- 
sistent with an independent calculation of the DDW 
modes starting from an effective low energy quantum ro- 
tor model, see Section lVll We note that the DDW modes 
found here may be overdamped, depending on the details 
of the system. This would be the case, for instance, if 
vbdw were smaller than the Fermi velocity. However, 
provided that the modes are underdamped, they can be 
detected in neutron scattering experiments. This is sur- 
prising at first, as local Zi invariance implies that equal- 
time correlators of the from (S™Sj) must decay expo- 
nentially at large distances |xj — Xj\. However, as (|2*7j) 
indicates, the correlator at different times need not de- 
cay exponentially. Also note that the inelastic scattering 
peak (|27|l has vanishing weight as k — > 0, consistent with 
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the absence of an elastic Bragg peak. Unfortunately, this 
feature makes DDW modes difficult to distinguish from 
the low energy collective modes of a quantum paramag- 
net, computed in Appendix iBl 



VI. QUANTUM ROTOR MODEL 

Another approach to study Goldstone modes is to in- 
troduce a quantum rotor model which captures the low 
energy properties of the system. This can be done by 
coarse- graining spins in plaquettes, as carried out in sec- 
tion HTD There we found that the low energy plaquette 
states at half-filling contain one singlet ground state, and 
one triplet state obtained by acting on the ground state 
by the staggered magnetization operator. This matches 
the low energy spectrum of a quantum rotor on the 
plaquette i, which can be written in terms of the compo- 
nent spins of the plaquette as 

L, = Sii + Si2 + Si3 + Si4, 
~ S,;i — Si2 + S^3 — Si4. 

Here, the total spin L acts as the canonical conjugate of 
n, satisfying the commutation relations, 



1 ' i _ 






i ' 71 j _ 






a 8 


= 0. 


(28) 



We introduce an effective Hamiltonian for these rotors, 
consistent with local Zi symmetry, 



(io) 



(29) 



The dynamic spin susceptibility can be obtained, in the 
long wave length limit k — > 0, by the linear response of 
L to an external field B, = Be- l ' krt - kT -', 



(30) 



The Heisenberg equations of motion for the rotors be- 
come 



The product of n operators then differs from the mean 
field result by a correction, 

n>f - %G tt/3 = 8 ljP al3 + X?f, 

which has been split into a local term p a @ , and a non-local 
term, A^ 3 = for i = j. Linearizing the equations of 
motion l|31(l with respect to these corrections, we obtain 

-iuoL a = F(k)e Q/37 G' 35 p 75 , 

-icjp af3 = ( t ^ s G aS + e a '' s G bd ){vL^ - B' 1 ) (32) 

where 



F(k) = 2g^(l-cos ft )- 



(33) 



Since we are interested in the linear response, we have 
excluded from Eq. Ij32(l terms quadratic in B. Note that, 
to this order in B, the non-local term drops out of 
the equations of motion. 

For the current case of interest, a uniaxial nematic, 
G a ° is of the form 



G = 



1 





n 2 - S 
n 2 + 2S 



Solving J2U for L yields 



L 1 



SF(k 2 ) 



lu 2 - isS 2 F(k) 



B 1 



(34) 



L n = 0, 



where _L is either of the directions perpendicular to n. 
Note that, unlike the case of antiferromagnetic order, x al3 
has no off-diagonal contributions. This is attributed to 
time reversal invariance. For long wave lengths k — > 0, 
F(k) — > gk 2 , and poles in the susceptibility indicate the 
presence of two degenerate DDW modes, with velocity 
«ddw = \J vS 2 g. This is consistent with the RPA result 
found above, Eq. (f27|l . 

Finally, we compute the dynamic spin susceptibility 
near k = (vt, it) by replacing the term if B by a source 
that couples directly to the staggered magnetization, 



HI 



(35) 



dhj 
dt 



jeNN(i) 



(hi ■ hj)(hi x hj) +B,xL, 



-CLi x n. ( - fij x Lj) - Bj x n. ( 



(31) 



In the absence of an external field B, the rotor correlator 
is assumed to be highly local, 

(hfh% - M^n 2 /3 + Q a0 ) = S t3 G af3 . 



Proceeding as above, we find the linear response 



-vG af} 



B p 

nr,2 % ' 



zgn 



(36) 



where z = 4 is the coordination number of the lattice. 
Note that, in the current approximation, the collective 
modes near (7r, it) are non-dispersing and gapped. The 
spin gap ui — y/ zgn 2 indicates that there is no long range 
Neel order in the system. 
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FIG. 9: (a) Phase diagram for eq. 1)3 7|l for the case J n = 0, 
corresponding to a Z2 lattice gauge theory with a U(l) field. 
The Ising transition as K is increased corresponds to a bind- 
ing of Q = 1/2 topological excitations into Q = 1 excitations, 
(b) Case Jg — 0, corresponding to a Z2 lattice gauge theory 
with an 50(3) field. The first order nematic to paramag- 
net transition can be split into two second order transitions, 
passing through a topologically ordered phase. 



VII. LATTICE GAUGE THEORY 
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Up to now, we have ignored the possibility of fraction- 
alization and the richer phase structure that it allows. 
As an example, by introducing a disclination core en- 
ergy, it is possible to split the phase transition for the 
onset of nematic order, which is strongly first order in 
Landau theory into two second order transitions [30| . 
For a theory with local Zi gauge redundancy, it is useful 
to discretize the magnetic degrees of freedom on a lattice 
and to introduce an auxiliary gauge field cry — ±1 living 
on the bonds of the lattice. Under a gauge transforma- 
tion at site i, the site variables pick up a minus sign, 
rii — > — rij and e l6i — > —e l9i , while simultaneously the 
bond variables surrounding i change sign, cry — ► — cry. 
The simplest gauge invariant action that can be written 
under these conditions is 

S = -J„ E njCTy-nj - J e ^ cry cos(6>i - dj) 

□ □ 

Unlike section llVl here we ignore the anisotropy in the 
tensor J"' 3 , which is weak and only affects the very low 
energy physics. Figures l§l and ITU1 display the phase di- 
agrams for extreme values of the couplings J n , Je and 
K. Figure |5| shows the phase diagrams for the cases 
J n = and Je — 0. These correspond to Z 2 lattice 
gauge theories with U(l) and 50(3) Higgs fields, re- 
spectively, which have been studied extensively in the 
literature 30, H3. 

Figure EH shows the situation when any of the cou- 
plings in (|37|l is infinite. In all of these cases, the aux- 
iliary field is completely ordered, and we can work 
in "minimal gauge", where ay — 1 for all bonds. This 
corresponds to a lack of topological defects, which are en- 
ergetically forbidden. Once a value of n (or 9) is chosen 
at a given lattice point, smoothness of the fields insures 



FIG. 10: Phase diagrams when any of the couplings in eq. 
13711 is infinite. In all of these cases, topological defects in the 
gauge field are energetically forbidden, and we can work in 
minimal gauge oy = 1. (a) Case J n = 00, insuring that at 
least nematic order is present. K is irrelevant in this situa- 
tion, as the field ay is frozen. As Jg is increased, the system 
undergoes an XY transition from nematic to SDW order. The 
situation is similar to the case Jg = 00, shown in (b), but now 
a Heisenberg transition is seen from CDW to SDW order. In 
(c), K = 00, leading to a dynamical decoupling of the 6 and n 
terms in the action 1371 . For small Jg and J n , a paramagnet 
with only topological order survives. 



that the Z2 redundancy in n (9) is removed everywhere. 
Thus, the universality class of these transitions is the 
same as that of the corresponding ungauged theories. 
For instance, in Fig. HOf c). starting with the topolog- 
ically ordered paramagnet at small J n and Jg, we can 
increase Je until the onset of CDW order through an 
XY phase transition. Similarly, by increasing J n we get 
the onset of nematic order through a Heisenberg transi- 
tion. In these regions of phase space, the simultaneous 
presence of both CDW and nematic orders imply SDW 
order. To see this, suppose that (n) m and {e l8 ) m are si- 
multaneously non-zero, where the subscript m indicates 
that we are working in minimal gauge. Then (e n) m is 
non-zero, but this quantity is in fact independent of the 
gauge used. 

Let us now consider the phase diagram when K = 0. 
In this case the auxiliary gauge field fluctuates strongly, 
and it is useful to sum l|37(l over configurations to get 

S eS = -f X>-^) 2 -f Ecos 2 ^-^) 
(ij) (y) 
- J / ^2n i -n j cos(9 i -9 J ) + ... (38) 
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FIG. 11: Five faces of the phase diagram of the action 137ft : 
Jg — 0, Jg — oo, J n = 0, J n — oo, and K — oo. We do 
not compute the phase diagram along the sixth face, K — 0, 
explicitly. However, there are indications from XY spins in 
2d that the only phases present there are paramagnet, CDW, 
SDW, and nematic0. However, introduction of a new term 
1391 can spit the SDW into a CDW+nematic phase. Note 
that, for large values of is', a topologically ordered phase be- 
comes available. 



In this particular model, J' = J n Je- An analysis of 
two dimensional algebraic order in a model of this sort 
(with XY spins) has been carried out in Ref. I34L There, 
the only phases present are paramagnet, CDW, nematic, 
and SDW, and there is a direct transition between para- 
magnetic and SDW phases. Thus, for the action Q37[l. 
although we only look at the boundary of the phase di- 
agram , in all the limits considered, the simultaneous 
appearance of nematic and CDW order implies SDW or- 
der. 

However, it is possible to have a state with nematic and 
CDW orders without simultaneously stabilizing SDW or- 
der. To see this, notice that J' in can be tuned by 
adding an extra gauge invariant term to the original ac- 
tion (j|7|) 



AS = — J e rii • rij cos(#2 — 9j). 



(39) 



In particular, when J' = 0, Eq. (|38J1 leads to two inde- 
pendent transitions for n and 0, yielding paramagnetic, 
CDW, nematic, and nematic+CDW phases. The lat- 
ter phase breaks SU(2) spin symmetry and translational 
symmetry, but leaves time reversal invariant. How can 
such a state arise from fluctuations of a spin stripe? In 
the stripe picture, an anti-phase domain wall serves a 
dual role, both as a region where charge accumulates 
and as a boundary between domains of opposite stag- 
gered magnetization. The effect of is to disentangle 
these two roles to obtain two separate objects, a charge 
line component and a magnetic anti-phase line compo- 
nent. Then, fluctuations of the anti-phase line compo- 
nent can restore time reversal symmetry without restor- 



ing translational invariance. Note that, as the extra cou- 
pling J e is decreased, an Ising phase transition between 
nematic+CDW and SDW phases begins to occur in the 
corner of the phase diagram Jg, J n — > oo until, for small 
enough J e , the nematic+CDW phase disappears entirely 
in favor of the SDW phase. 



VIII. SUMMARY 

In this paper we have considered the prospect of di- 
rect detection of a quantum spin nematic in strongly 
correlated electron systems, such as heavy-fermion com- 
pounds, the cuprates, and the organic superconductors. 
The spin nematic order parameter is a spin-2 operator, 
and it does not couple, to linear order, to many of the 
conventional probes, such as neutrons, photons, or nu- 
clear spins. However, we show that electrons moving in a 
nematic background have an anisotropic spin susceptibil- 
ity, which can be detected in polarized neutron scattering 
and Knight shift experiments. In addition, we discuss the 
possibility of observing the Goldstone modes associated 
with nematic ordering in inelastic neutron scattering ex- 
periments. 

In Section [HJ we defined the nematic order parame- 
ter in terms of the equal-time spin-correlation function, 
and argued for uniaxial nematic order in systems with 
collinear spin correlations. In Section IlIII we introduced 
wave functions to describe nematic order in a number 
of situations, including a stand-alone nematic, and ne- 
matic order coexisting either with spin-singlet or spin- 
triplet superconductivity. In order to do this, we fol- 
lowed a coarse-graining and low energy projection pro- 
cedure, as done in Refs. [z| and |2jJ We used these wave 
functions to constrain the short distance and low energy 
structure of the nematic order parameter. In Section Hvl 
we considered the finite temperature phase diagram of 
the nematic through a mean-field analysis of the GL 
free energy. Here, we included the effects of the spin- 
orbit interaction which, even when weak, ultimately fix 
the director of the nematic. We illustrated this princi- 
ple using spin-exchange constants measured in undoped 
La2Cu04. In Section we computed the anisotropic 
spin response of electrons in the nematic phase, and dis- 
cussed prospects for its observation. There, we also com- 
puted the spectrum of Goldstone modes of a quantum 
spin nematic. In Section lvTl these results were supported 
by the analysis of an effective quantum rotor model. Fi- 
nally, in Section IVIII we considered the possible pres- 
ence of topologically-ordered phases, as well as an exotic 
CDW+nematic phase, in a system of fractionalized elec- 
trons. 

Before concluding, we would like to mention a few ex- 
perimental systems where nematic order may be found 
(in addition to La a Sri_ 2; Cu04, which is discussed in the 
introduction [Hillll]). ^ n particular, heat capacity mea- 
surements in two-dimensional solid 3 He support the pres- 
ence of large many-spin interactions 35]. On the other 
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hand, numerical simulations of spins on a square lattice 
with ring exchange indicate a stable "p-nematic" phase 
In this case, frustration due to the ring exchange 
can destroy antiferromagnetic order in favor of a nematic 
phase. Another experimental system of interest is V2O3 
which, as temperature and pressure are varied, displays a 
Mott transition alon g a f irst order line ending in a second 
order critical point 37] . The lattice constants change 
discontinuously across the first order line, leading to a 
change in the ratio of potential to kinetic energies, and 
thus to a metal-insulator transition. The topology of 
the phase diagram is the same as that shown in the N~ 
phase surrounding point C in Fig.[fo] In fact, as discussed 
in section ITvl spin nematic ordering couples linearly to 
the anisotropy in the spin exchange. Thus, discontinuity 
in the spin nematic order parameter leads to a discontin- 
uous deformation of the lattice, which could explain the 
phase diagram of V2O3. 
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APPENDIX A: CORRELATIONS OF SPIN AND 
CHARGE 

While equation JSJ gives the simplest wave function 
that describes a state with both superconducting and 
nematic orders, it does not incorporate correlations be- 
tween charge and spin degrees of freedom. In the stripe 
picture, domain walls are not only regions where holes 
accumulate, they are also regions across which the stag- 
gered magnetization changes sign. One can modify © 
to include this effect, by relating Ui to the number of 
domain walls crossed: 



E IK 

{a 1 .<j 2 ...} = ±l k 



n 



(s + a im a il 



•<*J)|fi>< 



(Al) 



In this case, the sum over sites inside the delta function 
runs over a path on the plane joining site 1 to site k. This 
forces a change of sign in every time a new domain 
wall intervenes between those two sites. The path sum is 
only independent of path whenever the domain walls run 
continuously through the sample instead of ending up 



abruptly, and when an odd number of them do not inter- 
sect, as would happen, for instance, in a "T-junction". 
The functional g is chosen to destroy all configurations 
that violate these constraints, thus enforcing the integrity 
of the fluctuating domain walls. The simplest choice for 
g is a projection operator that gives equal weight to all 
allowed configurations. In this case, the wave function 
(IA1I) can be rewritten as 



|V>4 



EIK* 

V i<£D 



<rv ii m a & a J\il) i J] cSt|^, (A2) 



where T> denotes an allowed domain-wall configuration, 
and a-u.i is the sign of the staggered magnetization on 
site i in configuration T>. Equations (|A1I) and (|A2() are 



alternate descriptions of the nematic state of Zaanen and 
Nussinov 
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APPENDIX B: COLLECTIVE MODES OF A 
QUANTUM PARAMAGNET 

We would like to see if the director density waves 
(DDW) of a nematic state can be easily distinguished 
from the collective modes of a quantum paramagnet. 
These latter modes can be derived, in the RPA approxi- 
mation, by substituting —gQff — ► Jod a/3 gij in equations 
(|22|l and (|24|l . Here, r/ij is equal to 1 if i and j are 
nearest neighbors, and zero otherwise; and Jo > is the 
antiferromagnetic spin exchange. We find an isotropic 
spin response, 



XRP A ( k > W ) 



Xo(k,w) 



ra/3 



1 - J ry(k)xo(k,t 



(Bl) 



where ?y(k) = 2(cosfc l: + cosk y ). In the limit k — > 0, 



to — > 0, the free spin susceptibility tends to xo 
leading to the long wavelength result, 



fc 2 / W 2 , 



a(3 

Xrpa 



J 2 - 4Jnfc 2 



5 a(3_ 



Thus, there are three degenerate gapless paramagnon 
poles (with vanishing weight at k = 0). On the other 
hand, the poles of Eq. I|B1I) near k = (it, tt) have a spin 
gap A due to the lack of long-range Neel order, 



w k~(7T,7r) 



4J 2 (k~(7T^)) 2 +A2. 
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